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The qualities of electron refrigeration by means of tunnel junctions between superconducting and 
normal-metal electrodes are studied theoretically. A suitable approximation of the basic expression 
for the heat current across those tunnel junctions allows the investigation of several features of the 
device such as its optimal bias voltage, its maximal heat current, its optimal working point, and the 
maximally gained temperature reduction. Fortunately, the obtained results can be compared with 
those of a recent experiment. 
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I. INTRODUCTION 



The properties of the electron refrigeratoro are con- 
sidered in this paper. In conjunction with the SET- 
thermometera it forms a new devices for refrigerating 
and monitoring electrons at sub-Kelvin temperature. A 
future on-chip instrument can be used to provide isother- 
mal electrons to other experiments, e.g. in single electron- 
ics. The operation is similar to a household refrigerator in 
so far as the electron temperature, which decouples from 
the that of the phonons, can be chosen beforehand within 
a certain domain. Whereas the SET-thermometer, based 
on "orthodox theory" ,til is rather transparent, less effort 
is taken in case of the refrigerator .0 

The set-up of the experiment under considerationQ is 
shown schematically in Fig. [l| It consists of a pair of tun- 
nel junctions forming the refrigerator and another pair 
operating as thermometer. All junctions are connected 
to the island in the middle whose electrons are dealt with. 
The leads of the refrigerator are superconducting whereas 
the thermometer leads are made from a normal metal. 
The refrigerator double junction is biased by a voltage 
V close to A/e, where A is the superconducting gap of 
its leads. This allows only hot electrons (with an energy 
larger than the Fermi energy) to leave the island and, vice 
versa, only cold ones to enter it via the other junction. 
It is this dispersion which causes electron refrigeration at 
sub-Kelvin temperature, where the heating by phonons 
is weak. 

On the other hand, the thermometer represents an 
ultrasmall single-electron double tunnel junction. The 
thermal dependence of its Coulomb blockade enables the 
determination of the electron temperature. In general 
the thermometer is biased too. Since the zero-bias con- 
ductance is sufficiently and advantageous in order to de- 
termine the temperature, the bias is usually small and so 
is the thermometer current as well as the heating by this 
current. 

The tunneling regime of the refrigerator is considered 
in this paper, which is in contrast to Ref. ^. Due to 
the large junction resistance Andreev reflections are sup- 



pressed and therefore neglected. Following the argument 
from above, the heating by the thermometer current is 
also neglected. Furthermore, the refrigerating junctions 
are thought to be symmetric. This means that there 
is an equal voltage drop at each junction and the con- 
tinuity equatiou, for the electric current is fulfilled. In 
the experiment,!!! however, a slight asymmetry was found, 
which is neglected henceforth. Additionally, thermalized 
electrons are supposed, i.e. a sufficiently high temper- 
ature causes a sufficiently short electron-electron scat- 
tering time, which in turn causes thermalization of the 
electrons. As pointed out in Ref. || the electron tempera- 
ture should be above 100 mK to qualify this assumption. 
An opposite point of view,Q however, reveals interesting 
physics as well. 

Our approach starts from the classical formula for the 
heat current across a NIS junction, namelycl 



P (V) = ^J dsN(e)(e 



eV)[i N (e - eV) - f s {e)}, 



(1) 



where N(e) is the density of states on the supercon- 
ducting side, which is assumed to be the BCS den- 
sity of states, N(s) = |e|9(|e| - A)/Ve 2 - A 2 , and 
?N,s( £ ) = l/[ ex P( £ /^B^N,s) + 1] is the Fermi distribu- 
tion function for the normal- and superconducting side, 
respectively. R is the normal-state junction resistance. 
The integral (@) 

is solved in an analytic approximation in 
Sec. yl Later on conclusions of this result are discussed in 
Sec. jnj. It might arise the question, what the advantage 
of this approximative treatment is if a numerical solution 
is obtained readily. In our opinion those formulae yield a 
better insight into the behavior of the refrigerator and its 
qualities and therefore they are thoroughly worthwhile. 



II. CALCULATION 

In term of the BCS density of states N(e) Eq. [I] can 
be transformed into 
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P(V) 



J_ r d £ 

e 2 RJ A £ Ve 2 - A 2 
x[f N (e- e y) + f N ( e + e y)-2f s (e) 
_eV_ f°° e 
e 2 RJ A £ Ve 2 - A 2 

x[f N ( £ - e y)-f N (£ + ey)]. 



III. DISCUSSION 



(2) 



In order to simplify the calculation we^apply the following 
approximation to the Fermi function 



f( e )„e ( _ e ) + !^f)e->M 



(3) 



with 7 = l/(2fc B Tln2) and sgn(e) = e/\e\, sgn(O) = 0. 
The value of 7 is adjusted to meet 



J def(e) = - jf de< 



In the limit T — > the approximation is exact. Fur- 
thermore, it is right at s = 0. For large |e| ^> k^T the 
difference to the exact formula is negligible, but for an 
intermediate energy, e ~ k&T it is not. Hence, we expect 
this approximation to work out at low temperature, i.e. 
at kftT <C A here,-. This condition is obeyed fairly good 
in the experiments 

The main advantage of (0), however, is the possibility 
to reduce the integrals of (|2j) to Laplace transformations. 
In detail we find with 

f N (e -eV) + f N (e + eV) - 2f s (e) « e" 7N£ cosh( 7N e V) 

-e- 7s£ , 

f N (e - e V) - f N (e + e 7) w e" 7NE sinh(7 N e 1/) 



ant 



dt/(t)e 



-pt 



e9(e - A) 

Ve 2 - A 2 
e 9(e - A) 
Ve 2 - A 2 



( 7 ) = AK 1 ( 7 A), 



( 7 ) = -A— Ki( 7 A) 

= ^-[K ( 7 A) + K 2 ( 7 A)], 



where K l/ (z) is the modified Bcssel function of ^-th or- 
der, 



P(V) 



A 2 

^-^{[K (7nA) + K 2 ( 7N A)] cosh( 7N e V) 



-K (7 S A)-K 2 ( 7s A)} 
Ae Vsinh(7Ne V) 



e 2 R 



Ki(tnA). 



(4) 



Eq. [| is the main result of this paper. We are going to 
discuss it in the next section. 



Following Ref. ^| we consider the limit 7n,sA ^> 1, 
e V = A first, which describes the experimental situation 
to a large,-extent. In this case the use of the asymptotic 

M 



expansions 



K v (z) 



simplifies (0) to 



4^ 2 -l 



P(A/e 



8z 



(7rln2) 3 / 2 A 2 fk B T N \ 
2vr e 2 R\ A J 

e 2 R V A 



3/2 



(5) 



Eq. H recovers the temperature dependence P(A/e) oc 

3/2 13 
T N of the exact result .□ The numerical prefactor, how- 
ever, is found to be slightly larger (0.48 in Ref. ^). The 
agreement approves of the introduced approximation (|^) 
and indicates that its accuracy is in the few-percent 
range. 

In a next step (Q) can be used to determine the optimal 
value of the bias voltage V . From the formal derivative 
of (Q) it is found that this value V opt should obey the 
equation 

T/ fW T/ \ a dln[K 1 ( 7N A)] 
7Ne Vopt coth(7 N e Vopt) = -7nA — 1. 



d( 7N A) 



(6) 



We like to stress that this result is independent of the 
lead temperature Tg. In the limit 7nA 3> 1 and e V ~ A 
Eq. simplifies again considerably to 



eV 



pt 



A - fc B T N In 2 w A - 0.69 /c B T N 



(7) 



In Fig. ||we plot both the numerical solution of (^J) and its 
approximation ^7§. Up to k^T^ ~ 0.3 A the approxima- 
tion works nicely thus covering the relevant temperature 
range (see next paragraph). It is difficult to observe this 
shift in the experimental datajj since it corresponds to 
maximal 28 yiV only. 

In Fig. |^ the heat current (0) using the optimal bias 
voltage V"opt of @ is shown in comparison with both ex- 
perimental data and the simpler approximation (^) . The 
maximum of the heat current is found at &bTn < 0.23A 
always. Therefore (Q) is a practical approximation under 
relevant conditions as outlined above. Whereas Eq. |^ 
yields a useful rule of thumb, the numerical solution is 
rather convincing. Both the optimal temperature and 
the maximal heat current are determined satisfactory for 
higher temperature. For lower temperature, however, we 
observe an increasing deviation. This deviation can not 
be due to our approximation (j^), which improves with 
decreasing temperature. According to Ref. || also the 
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thermometer as such is not responsible for those discrep- 
ancies. We attribute the deviation to rising nonequilib- 
rium effects,! 2 ! i.e. a non-Fermi like energy distribution of 
the electrons on the island. In Rcf. || a value of 100 mK 
is given as a lower limit of the equilibrium description. 
The exact value, however, depends on the specific sam- 
ple layout. It is reasonable to expect an increased limit 
for the experiment under consideration since the temper- 
ature reduction is considerably larger and consequently 
the nonequilibrium effect becomes more pronounced. 

Finally, we are going to discuss an estimate of the lower 
limit of the electron temperature Tn based on (jij). It is 
assumed that the heating of the electrons cp, the island 
is only due to phonon couplingu and followst^! 



Pei-ph = SQ(T| - T N ) 



(8) 



with the island volume f2 and a material dependent pref- 
actor E. Eq. |s| implies the assumption that the electrons 
in the leads are in thermal equilibrium with the phonons 
and the phonon temperature is constant throughout the 
whole device. For the stationary situation at optimal 
bias, P(V op t) — Pei-ph^ follows a relation between Tn 
and Tg in terms of Eqs. [|, |[ and 0. Simpler is the treat- 
ment using (||) and (0) instead of (Q) and (^) , respectively, 
which results in the expression 



A 



fc B T N 



,4 



3/2 



(9) 



with A = (7rln2) 3 / 2 /(27r)A 2 /(e 2 i?Erj)(fc B /A) 5 . For the 
maximum temperature reduction, Ts — Tn, (|9|) yields 



opt 



A 



feT^opt 



0.3A 



kB?N,opt 



1/2 



(10) 



for the optimal values of Tn,s- Using a further approxi- 
mation, T Ni0p t <§; Ts,o P t, we find from @ and (|i"o| ) 



TN,o P t«0.3 10 / 7 A 2 / 7 «0.179A 2 / 7 , 
Ts,o P t ~ 0.3 3/7 A 2/7 w 0.597^ 2/7 . 



(11) 



The last approximation reveals an interesting uniform 



relation T 



N,opt/T3,opt 



0.3 for the optimal working 
point of the refrigerator, independent of all device pa- 
rameters. Indeed, the optimal values of the experiment, 
Tk.opt ~ 100 mK and Is, opt ~ 300 mK, almost meet this 
universal ratio. Furthermore, since the figure of the ratio 
assembles from the general exponents 3/2 and 5 only, it 
is expected to be right beyond the range of our approxi- 
mation. 

In Fig. [|thc numerical solution of P(V op t) = Pi-ph in 
terms of Eqs. ^, ^, and |^ is compared with the approxi- 
mation (^) and the experimental data. Similarly to the 
discussion above we attribute the deviations at low tem- 
perature to the manifestation of nonequilibrium. Other- 
wise we notice fair agreement of the numerical solution 



and its approximation as long as the temperature is not 
to high. In comparison with the experiment we compute 
?kopt « 75 mK and T S:0p t « 250 mK from flu]). This 
values are somewhat lower than the measured data (see 
last paragraph). 



IV. CONCLUSION 



In this paper electron refrigeration by means of tunnel 
junctions is studied theoretically. An approximation of 
the original expression of the heat current across those 
junctions enables a closed analytic formula for this quan- 
tity. Based on this formula several features of the refrig- 
erator are discussed. The optimal bias voltage is deter- 
mined and the proportionality of the corresponding heat 
current to T^ 2 is confirmed. The maximal difference 
of the electron temperature in the leads, Ts, and the is- 
land, Tm, is discussed in stationary regime. Analytic ap- 
proximations to the optimal values of these temperatures 
are derived and an universal ratio TN, pt/Ts, pt = 0.3 is 
found. All results are discussed in comparison with a 
recent experiment. The investigation of nonequilibrium 
effects, that occur at very low temperature, remains an 
interesting open problem. 
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FIG. 1. Combination of refrigerator and thermome- 
ter, both based upon tunnel junctions. The hori- 
zontally shown refrigerator consists of a superconduc- 
tor-insulator-normal-metal (SIN) junction to the left and a 
NIS-junction to the right biased by a voltage V. The verti- 
cally displayed thermometer tests the electron temperature of 
the island in the middle by means of two NIN-junctions. 
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FIG. 3. The maximal heat current P(V op t) as a function 
of the electron temperature Tn. The solid lines correspond 
to different ratios Ts/Tk of the electron temperature in the 
leads and the island. For the uppermost curve this ratio is 
1, whereas it corresponds to the values of the experiment for 
the other curves. The dashed line shows Eq. ^, and the dots 
(•) correspond to the values of the experiment. 
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FIG. 2. Optimal bias voltage V op t in dependence on the 
electron temperature on the island Tn. The solid curve dis- 
plays the numerical solution of (Q) and the dashed line its 
approximation (M). 



FIG. 4. Temperature reduction Ts — Tn versus Tn. The 
solid curve corresponds to the numerical procedure described 
in the text, whereas the dashed line is a plot of (^). The 
symbols (•) are again the data of the experiment. 
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